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ABSTRACT 

Solar activity has significantly changed over the last two Schwabe cycles. After a long and deep minimum at the end of Cycle 23, the 
weaker activity of Cycle 24 contrasts with the previous cycles. In this work, the response of the solar acoustic oscillations to solar 
activity is used in order to provide insights on the structural and magnetic changes in the sub-surface layers of the Sun during this 
on-going unusual period of low activity. We analyze 18 years of continuous observations of the solar acoustic oscillations collected by 
the Sun-as-a-star GOLF instrument onboard the SoHO spacecraft. From the fitted mode frequencies, the temporal variability of the 
frequency shifts of the radial, dipolar, and quadrupolar modes are studied for different frequency ranges which are sensitive to different 
layers in the solar sub-surface interior. The low-frequency modes show nearly unchanged frequency shifts between Cycles 23 and 24, 
with a time evolving signature of the quasi-biennial oscillation, which is particularly visible for the quadrupole component revealing 
the presence of a complex magnetic structure. The modes at higher frequencies show frequency shifts 30% smaller during Cycle 24, 
which is in agreement with the decrease observed in the surface activity between Cycles 23 and 24. The analysis of 18 years of GOLF 
oscillations indicates that the structural and magnetic changes responsible for the frequency shifts remained comparable between 
Cycle 23 and Cycle 24 in the deeper sub-surface layers below 1400 km as revealed by the low-frequency modes. The frequency shifts 
of the higher-frequency modes, sensitive to shallower regions, show that Cycle 24 is magnetically weaker in the upper layers of Sun. 
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1. Introduction 


The unexpected long activity minimum of the solar Cycle 23 
in 2007-2008 puzzled over the complexity and the variability 
of the emerging solar magnetic fields generated in the convec¬ 
tive zone. It revealed even more the difficulty to perform reliable 
prediction of the 11-year solar cycle. The current Cycle 24 is 
showing one of the weakest level of activity since a century, with 
a decreased s urface activity by about 30% compared to Cycle 
23. Moreover. iLivingston et al.l (1201 2l) observed a continue drop 
in sunspot magnetic field strength since the last 30 years argu¬ 
ing that sunspots could disappear from the Sun’s surface during 
Cycle 25 as the 1500 G threshold to form dark sunsp ots would 
be reached. Deeper inside the Sun. lHowe et al.l (1201 3l) have used 
helioseismic data to follow the zonal flow pattern, called tor¬ 
sional oscillation, during Cycles 23 and 24. They showed in par¬ 
ticular that the high-latitudinal pattern, one of the signatures of 
a new Cycle, is actually absent since the beginning of Cycle 24, 
probably due to the lower polar magnetic field strength. 


The solar acoustic (p) oscillation frequencies were revealed 
to respond to changes in the surface activity, and to vary in 
correlation with solar proxies such the sunspot number or the 
10.7-cm radio emission. The temporal variations of the so¬ 
lar P 2 mode_freguencies_w^ solar activity were first reported 
^ Woodard & NovesI d 19851). shortly co nfirmed by iFossat et"^ 
(Il987h and lLibbrecht & Woodardl(ll990l) . They were quickly at¬ 
tributed to changes taking place in the sub-surface layers. It is 
also known for more than a decade that the extraction of the 
sound-speed profile in the core must be done carefully. It is why 
this information was obtained in using only low-degree modes 


with frequencies smaller than 2.4 mHz (iTurck-Chieze et aTL 
IMTtiTurck-Chieze & Lo^l20H . 

The solar frequencies are today identified as a unique proxy 
of solar activity, as they might reveal inferences on sub-surface 
changes with solar activity not visible at the surface. They could 
also contribute to build varying topologies of magnetic field 
below the surface. As longer and higher-quality helioseismic 
observations became available, the temporal variability of the 
p-mode frequencies was studied in more details. These fre¬ 
quency shifts were actually observed to be frequency depen¬ 
dent, the shifts being larger at higher frequencies, and to be 
angular-degree (1) dependent, i.e., mode-inertia dependent, al- 
tho ugh the / dependence is small for low-degree modes (see 
e.g. Gellv et al.L l2 002t Howe et al.l l2002t iJiminez-Reves et al.L 


12004 : Salabert et al. r i2004l) . The shifts can be then interpreted 

to arise from structural chang es in t h e sub -surface layers, such 
as changes in tem perature (iKuh nl 19881) or in the size of 
the acoustic cavity (iDziembowski & Goodel l2005l) . Moreover, 
iHowe et ^ (l2002h showed that the temporal and latitudinal dis¬ 
tribution of the frequency shifts is correlated with the spatial dis¬ 
tribution of the surface magnetic field. Thus, the frequency shifts 
cannot be purely explained by structural changes and the deter¬ 
mination of the dipolar and the toroidal magnetic fields of the 
sub-surface layers appear as one of the important outputs from 
these observables (iBaldner et al.Ll2()()^ . 

However, while the frequency shifts varied closely to so¬ 
lar surface activity with very high levels of correlation over 
the last three complete solar Cy cles (21, 22, and 23 ), span¬ 
ning thirty years of solar activity (IChaplin et al.L l2007h . signif¬ 
icant differences were observed between acoustic frequencies 
and solar activity during the unusual minimum of Cycle 23 as 
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reported in Salabert et alJ (l2009h . Broomhd _et^ (l2009t) . and 
iTripathv et al.l (l201(ll) . Moreover. iBasu et alJ ( 2012h showed that 
the structure of the solar sub-surface layers was very different be¬ 
tween Cycle 22 and Cycle 23, with deeper changes during Cycle 
22 compared to Cycle 23. 

In this study we investigate in more details the differences in 
the variability of the radial, dipolar, and quadrupolar low-degree 
modes. We show how the temporal evolution of the magnetic 
configuration imprints differently these modes during Cycle 23 
and the weak on-going Cycle 24 by analyzing the low-degree 
oscillation frequencies extracted from the observations collected 
by the space-based, Sun-as-a -star Global Oscillat ions at Low 
Frequency instrument (GOLF: [Gabriel et al.Lll995h onboard the 
Solar and Heliospheric Observatory (SoHoT IPomingo et al.L 
Il995h spacecraft. SoHO was launched in December 2, 1995, 
and since then is collecting continuous data of the Sun from a 
stable orbit around the LI Lagrangian point. With more than 
18 years of space observations available today, free from any 
perturbations related to the Earth atmosphere as it occurs from 
ground-based observations, the GOLF/SoHO observations 
provide an exquisite and unique dataset for global helioseismic 
studies of the low-degree oscillation modes. In Section 2, 
we describe the data used in this paper and how the seismic 
parameters were extracted. The frequency shifts thus estimated 
are presented and discussed in Section 3. In Section 4, we 
provide possible explanations and interpretations of the related 
structural and magnetic changes of the sub-surface layers 
during Cycle 24. The main results are summarized in Section 5. 
Finally, the tables containing the oscillation frequencies of the 
GOLF low-degree modes extracted from the 365-day spectra 
analyzed in this study are made electronically available as 
described in Appendix A. In Appendix B, 18 years of frequency 
shifts at individual angular degree I are presented in more details. 


2. Data and analysis 

2.1. Data sets 

We analyzed 18 years of continuous observations collected by 
the GOLF instrument, a resonant scattering spectrophotometer 
measuring the Doppler wavelength shift - integrated over the 
solar surface - in the D1 and D2 Fraunhofer sodium lines at 
589.6 and 589.0 nm respectively. It has been designed to mea¬ 
sure the Doppler shift by switching from one side of the wing 
to the other side, an additional coil allowing two near point mea¬ 
surements along_eachsideofthe wings in the more stable part of 
the lines dGabriel et al.L Il9^ . However, the GOLF instrument 
has not been used in its nominal configuration due to a malfunc¬ 
tion in its polari zation switching m echanism identified shortly 
after the launch (iGarcfa et al.L l2004li . As a consequence, obser¬ 
vations from only one side of the sodium doublet are collected, 
from which a proxy for the Doppler velocity signal is formed. 
Moreover, observations have been collected from each side of 
the doublet as follows: in the blue-wing configuration from 1996 
April 11 until 1998 June 26, and later on from 2002 November 
19 until now (the so-called blue periods)', in the red-wing con¬ 
figuration of the sodium doublet between 1998 October 30 until 
2002 November 18 (the so-called red period). These two con- 
figurations imply d ifferent spatial weightings of the solar disk 
(iGarcia et al.Lll99^ . 

In addition, two extended gaps are present in the analyzed 
time series because of the temporary loss of the SoHO space¬ 
craft. The first gap of about 100 days happened during the sum¬ 


mer of 1998 after a bad maneuver of the SoHO rotation space¬ 
craft. The second one of a period around 1 month in January 
1999 occurred while new software was being uploaded to the 
spacecraft due to a failure in the SoHO gyroscopes. 

The observational blue-wing configuration of the GOLF in¬ 
strument has not been modified since November 2002, except 
for readjusting the photomultiplier high voltages in order to 
compensate its aging. It corresponds today to almost 12 years 
of continuous measurements, with a duty cycle close to 100%. 
The response function of the GOLF sodium Fraunhofer line 
is different between t he blue- and the red-wing configurations 
(iLefebvre et al.L 1^081) . with averaged hei ghts in the solar atmo¬ 
spher e of 322 km and 480 km respectively Oimenez-Reves et all 
l2007h . Then, the measurements obtained in the red-wing config¬ 
uration originate higher up in the solar atmosphere with larger 
contributions from the chromosphere. 

A total of 6538 days of GOLF observations with a duty 
cycle of 96.0% covering 18 years between April 11, 1996 
and March 5, 2014 were analyz ed. The GOLF velo city time 
series were obtained following iGarcfa et al.l (l2005l) and the 
a mplitudes of the blue - and r ed-wing signals normalized as 
in iJimenez-Reves et ^ (l2003l) . To obtain relations between 
mode frequencies and solar activity, this dataset was split into 
70 contiguous 365-day sub series, with a four-time overlap of 
91.25 days. 


2.2. Extraction of the p-mode parameters 

The power spectrum of each time series was fitted to estimate the 
mode parameters of the I = 0, 1,2, and 3 modes as described in 
[Salabert et al.l (l2007l) . Each mode component of radial order n, 
angular degree I, and azimuthal or der m was parameterized wit h 
an asymmetric Lorentzian profile (iNigam & Kosovichevlll99^ . 
as: 


dln.i.mix) — Hi,i 


(1 -I- bnjXnj)^ + blj 

1+X^, 

n,l 


( 1 ) 


where x„j - 2(v - v„j)IT„j, and v„j, F„_/, and H„j represent the 
mode frequency, linewidth, and height of the spectral density, 
respectively. The peak asymmetry is described by the parameter 
b„ i. Because of their close proximity in frequency, modes are 
fitted in pairs (i.e., I — 2 with 0, and I - 3 with 1). While 
each mode parameter within a pair of modes are free, the peak 
asymmetry is set to be the same within pairs of modes. When 
being present inside the fitting window, the / = 4 and 5 modes 
were included in the fitted model. An additive value B is added 
to the fitted profile representing a constant background noise in 
the fitted window. Since SoHO observes the Sun equatorwards, 
only the I -i- \m\ even components are visible in Sun-as-a-star 
observations from the GOLF instrument. The amplitude ratios 
between the I = 0,1,2, and 3 modes and the m-height r a tios of 
the I - 2 and 3 multiplets calculated in ISalabert et al.l (1201 ih 
for the GOLF measurements were used. Finally, the mode pa¬ 
rameters were extracted by maximizing the likelihood function, 
the power spectrum statistics being by a with two degrees 
of freedom distribution. The natural logarithms of the mode 
height, linewidth, and background noise were varied resulting in 
normal distributions. The formal uncertainties in each parameter 
were then derived from the inverse Hessian matrix. 
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1995 2001 2006 2012 
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Fig. 1. Temporal variations of the frequency shifts in //Hz averaged over the modes 1 — 0,1, and 2, (Av„_;=o,i, 2 ), and calculated for 
four different frequency ranges (black dots). From left to right, the frequency ranges are the following: a) 1800 //Hz < v < 3790//Hz; 
b) 1800 //Hz < V < 2450 //Hz; c) 2450 //Hz < v < 3110 //Hz; and d) 3110 //Hz < v < 3790 //Hz. The 10.7-cm radio flux, F 10 . 7 , 
averaged over the same 365-day timespan and scaled to match the rising phase and maximum of Cycle 23 is shown as a proxy of 
the solar surface activity (solid line). 
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Fig. 2. Same as Fig. [T]but after removing the QBO’s signature from the temporal variations of the frequency shifts, (Av„/=o,i, 2 ) 
(black dots). The 10.7-cm radio flux, Fiqj, was smoothed in the same way and scaled to match the rising phase and maximum of 
Cycle 23 (solid line). 


3. Results 

The temporal variations of the mode frequencies, {Av„j), were 
defined as the differences between the frequencies observed at 
different dates and reference values of the corresponding modes. 
The set of reference frequencies was taken as the average over 
the years 1996-1997, corresponding to the minimum of the 
activity Cycle 22. The formal uncertainties resulting from the 
peak-fitting analysis were used as weights in the average com¬ 
putation. The frequency shifts thus obtained were then weigthed 
averaged over 15 consecutive radial orders between 1800 //Hz 
and 3790 //Hz. The lower frequency was defined by the lowest 
modes which are accurately fitted in all the analyzed 365-day 
spectra, i.e., n = 12 for / = 0 and 1, and n = 11 for / = 2. Note 
that we did not use the / = 3 mode in the following analysis 
because of its lower signal-to-noise ratio. In addition, mean val¬ 
ues of daily measurements of the 10.7-cm radio flux, F 10 . 7 , were 
used as a proxy of the solar surface activitjQ. 


3.1. At different depths below the surface 


As demonstrated 


Howe et al. 

(|l999 

iJimenez-Reves et al. 

Salabert et alJ 

(|200^ 

]) for both low- and 


letakl (Il998h . 
(I2001h . and 
and medium-angular 


* The 10.7-cm radio flux data are available 
from the National Geophysical Data Center at 
http://www.ngdc.noaa.gov/stp/solar/solardataservices.html 


degrees, the frequency shifts associated to the solar Schwabe 
cycle present a strong frequency dependence, going from small 
shifts at low frequency to much larger shifts at high frequency. 
As modes with different frequencies are sensitive to different 
layers in the sub-surface outer parts of the Sun, this information 
can be used to examine differences between solar Cycles 23 and 
24 at different depths. 

We thus defined three distinct frequency ranges containing 
each five consecutive radial orders covering different reflect¬ 
ing points in the solar atmosphere. The three frequency ranges 
were defined as follows: (1) the low-frequency range where 
1800 //Hz < V < 2450 //Hz; (2) the medium-frequency range 
where 2450 //Hz < v < 3110 //Hz; and (3) the high-frequency 
range where 3110 //Hz < v < 3790 //Hz. These three frequency 
ranges allow us to study the response of the p-mode oscillations 
to solar activity at different depths below the surface, and thus 
to infer localizations of the layers where the magne t ic field s af¬ 
fect the oscillation frequencies. Indeed, iBasu et akl (l20I2h cal¬ 
culated that the averaged kernels of the mid- and the high- 
frequency ranges have their largest sensitivities at 0.9981 Rq 
(i.e., 1300 km) and 0.9989 Rq (i.e., 760 km) respectively, while 
the low-frequency modes peak deeper in the solar interior at 
0.9965 Rq (i.e., 2400 km). Nevertheless, while the mid- and the 
high-frequency modes correspond to a thin layer below the Sun’s 
surface, the averaged kernel of the low-frequency modes covers 
a much wider region ranging from 1400 km down to 4000 km 
deep. Finally, four frequency ranges were thus used to calculate 
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Table 1. Frequency shifts per unit of change in the 10.7-cm solar radio flux and linear correlations extracted from the analysis 
of independent 365-day GOLF spectra for different phases of the solar Cycles 23 and 24. The values were obtained from the 15 
consecutive radial orders between 1800 //Hz < v < 3790 //Hz, and for 5 consecutive radial orders in the low-frequency range 
1800 //Hz < V < 2450 //Hz, 5 consecutive radial orders in the mid-frequency range 2450 //Hz < v < 3110 //Hz, and 5 consecutive 
radial orders in the high-frequency range 3110 //Hz < v < 3790 //Hz. 



All data 

Rising phase Cycle 23 

Declining phase Cycle 23 

Rising phase Cycle 24 


(Apr. 1996 / Apr. 2014) 

(Apr. 1996/Oct. 2001) 

(Oct. 2001/Jan. 2009) 

(Jan. 2009 / Apr. 2014) 

1 

Gradient r,, 

Gradient" 

Gradient" r,, 

Gradient" r,, 


1800//Hz< v< 3790//Hz 


(0,1,2) 

1.7 ± 0.1 

0.95 

1.4 ±0.1 

0.98 


1.6 ±0.1 

0.95 

2.1 ±0.1 

0.95 

0 

1.7 ± 0.1 

0.95 

1.2 ±0.1 

0.97 


1.7 ± 0.1 

0.97 

2.1 ±0.2 

0.99 

1 

2.1 ± 0.1 

0.95 

1.7 ±0.1 

0.97 


1.7 ± 0.1 

0.95 

2.3 ± 0.2 

0.86 

2 

1.4 ±0.1 

0.82 

1.4 ±0.1 

0.94 


1.3 ±0.1 

0.81 

2.0 ± 0.2 

0.75 





1800 //Hz 

< 

V < 2450 //Hz 




(0,1,2) 

0.4 ±0.1 

0.49 

0.4 ±0.1 

0.77 


0.3 ±0.1 

0.50 

0.8 ± 0.3 

0.56 

0 

0.5 ± 0.2 

0.38 

0.6 ±0.1 

0.74 


0.5 ± 0.2 

0.58 

0.5 ± 0.3 

0.33 

1 

0.4 ± 0.2 

0.46 

0.3 ±0.1 

0.59 


0.3 ±0.1 

0.26 

0.9 ± 0.3 

0.48 

2 

0.4 ±0.1 

0.33 

0.4 ±0.1 

0.54 


0.1 ±0.1 

0.13 

1.0 ±0.3 

0.40 





2450 //Hz 

< 

y < 3110 //Hz 




(0,1,2) 

2.1 ±0.1 

0.97 

2.0 ±0.1 

0.96 


2.3 ±0.1 

0.95 

2.4 ± 0.2 

0.96 

0 

1.5 ±0.2 

0.92 

1.1 ±0.1 

0.95 


1.9 ±0.1 

0.95 

2.0 ± 0.2 

0.94 

1 

2.5 ± 0.2 

0.97 

2.7 ±0.1 

0.98 


2.5 ±0.1 

0.97 

2.4 ± 0.3 

0.87 

2 

2.3 ± 0.2 

0.89 

2.3 ± 0.2 

0.96 


2.8 ± 0.2 

0.89 

3.2 ±0.3 

0.91 





3110//Hz 

< 

V < 3790 //Hz 




(0,1,2) 

4.6 ± 0.2 

0.98 

4.6 ±0.1 

0.99 


4.6 ± 0.1 

0.95 

5.2 ±0.3 

0.99 

0 

4.1 ±0.3 

0.96 

3.8 ±0.2 

0.97 


4.2 ± 0.2 

0.96 

5.5 ± 0.4 

0.96 

1 

5.2 ±0.3 

0.97 

5.3 ±0.2 

0.99 


5.2 ± 0.3 

0.93 

5.7 ± 0.4 

0.96 

2 

4.7 ± 0.3 

0.94 

4.9 ± 0.3 

0.95 


4.6 ± 0.3 

0.95 

4.1 ±0.5 

0.93 


“ Gradient against the 10.7-cm radio flux in units of nHz RF ‘ (with 1 RF = 10 J s ‘ m ^ Hz *). 


the frequency shifts, when we include the one between 1800 //Hz 
and 3790 //Hz. 

Figure [T] shows the temporal variations of the frequency 
shifts, (Av„y=o,i, 2 ), averaged over the modes / = 0,1, and 2 and 
calculated over the four defined frequency ranges. The corre¬ 
sponding 10.7-cm radio flux averaged over the same 365-day sub 
series is also shown for comparison. It was scaled to match the 
rising phase and the maximum of the solar Cycle 23, by perform¬ 
ing a linear least-square regression with the frequency shifts. In 
Table[Tl we present the frequency shifts, (Av„j=o,i, 2 ), calculated 
per unit of change in the 10.7-cm radio flux obtained by mini¬ 
mizing weighted linear regressions between these two quantities, 
as well as the associated linear correlation coefficients, derived 
for the different frequency ranges. Note also that the values given 
in Table [1] correspond to one set of independent spectra (one 
over four points) and that they are consistent within Icr with any 
other of the chosen sets. As already widely established (see e.g. 
the references given in Section 1), the temporal variabilities of 
the p-mode oscillation frequencies extracted from 18 years of 
GOLF data are observed to be closely correlated with surface 
activity following the 11-year periodicity of the solar cycle. The 
correlation coefficients are over 0.95 for the mid- and the high- 
frequency ranges, but the low-frequency modes show a smaller 
correlation with surface activity of about 0.5. While the higher 
frequency modes show the largest variations along the solar cy¬ 
cle with frequency shifts being more than 10 times larger than 
the ones at low frequency, the low-frequency p modes still show 
significant variations between periods of low and high surface 
activity (Fig.[TJ. Moreover, the minimum of Cycle 23 in the fre¬ 
quency shifts of the mid- and high-frequency modes is observed 


to be lower than the minimum of Cycle 22 in May-June 1996, 
which is in agreement with the weaker surface activity since the 
long and deep minimum of Cycle 23 compared to Cycle 22. The 
quasi-biennial oscillation (QBO) (iFletcher et al.L 1201 Oh with an 
amplitude modulated by the 11-year solar cycle is also visible 
on Fig.[T]in all the frequency ranges. 

3.2. At different periods of the cycles 

Significant differences in the relationship between acoustic os¬ 
cillation frequencies and activity proxies were observed dur¬ 
ing the unusual extende d minimum of Cycle 23 compared to 
other phases of activ i ty (ISalabert et al.L l2009t iBroomhall et al.L 
120091: l^ipathv et ^ l2010h . For these reasons, we examined 
three separate periods: a) the rising phase of Cycle 23 going 
from April 1996 to October 2001 with a mean 10.7-cm radio 
flux of 129 RF3; b) the declining phase of Cycle 23 going from 
October 2001 to January 2009 with a mean 10.7-cm radio flux of 
114 RF; and c) the rising phase of Cycle 24 going from January 
2009 to April 2014 with a mean 10.7-cm radio flux of 99 RF. 
The starting and ending dates of each of the three activity phases 
were defined using the 10.7-cm radio flux properl y smoothed to 
remov e the signature of the QBO as described in iFletcher et al.l 
(l2010h . Table [1] gives the gradients and the correlation coeffi¬ 
cients between frequency shifts and radio flux during the three 
phases of activity for the four frequency ranges. It shows that 
the frequency shifts averaged over 1800 //Hz < v < 3790 //Hz 
are observed to be about 50% larger per unit of change in radio 


^ The 10.7-cm radio flux has for units 1 RF = 10 J s ' m ^ Hz '. 
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flux during the rise of Cycle 24 than during the rise of Cycle 
23, whereas the radio emission decreased by about 30% over the 
same period of time. In addition, during the declining phase of 
Cycle 23, the gradients are observed to be about 15% larger than 
during the rising phase of Cycle 23. Moreover, the discrepancy 
in the relationship between frequencies and radio flux is more 
important for the low-frequency modes, i.e., for the modes go¬ 
ing into deeper layers of the Sun, than for the higher-frequency 
modes confined in shallower layers, and it is accentuated during 
the rise of Cycle 24 (Table [TJ. The correlation coefficients with 
solar activity indicate as well a decrease of almost 40% between 
Cycle 23 and Cycle 24 for the low-frequency modes, while for 
the mid- and the high-frequency modes they remained similar. 

Figure |2] shows the frequency shifts, (Av„y=o,i, 2 ), averaged 
over the modes / = 0,1, and 2 once the QBO’s signature was re- 
m oved by applying a pr oper smoothing of 2.5 years as prescribed 
in iFletcher et alJ (1201 Ol) for the four frequency ranges previously 
defined. The solid line on each panel corresponds to the scaled 
10.7-cm radio flux, smoothed in the same manner. Note that in 
order to avoid misinterpretation due to possible border effects 
introduced by the smoothing, we do not consider the extreme 
points corresponding to the first three and the last three points 
of the smoothed data. The larger frequency shifts compared to 
the surface activity proxy is clearly seen on Fig. |2] during the 
rise of Cycle 24 from 2008. The low-frequency modes show 
the largest deviation with the 10.7-cm radio flux, with much 
larger frequency shifts than would be expected based on the sur¬ 
face activity proxy. The declining phase associated to Cycle 23 
is much steeper as well than the activity proxy, unlike what is 
observed for the mid- and high-frequency p modes. Moreover, 
while the mid- and high-frequency ranges show temporal varia¬ 
tions cleaned from the QBO, a longer period seems to appear 
in the low-frequency modes since 2006. Although that needs 
further observations in order to be confirmed, it could indicate 
changes related to the QBO deeper than 1400 km in the Sun’s 
sub-surface layers. See Appendix B for additional figures based 
on the analysis of the frequency shifts at individual angular de¬ 
gree. 

The left panel of Fig. 0 shows the frequency shifts, 
(Av„/=o,i, 2 ), averaged over the 15 radial orders between 
1800 yuHz and 3790 /rHz as a function of the associated 10.7-cm 
radio flux for the three phases of solar activity defined above, 
both q uantities being properly smoothed as in iFletcher et alJ 
(120101) . A clear magnetic hysteresis appears between oscilla¬ 
tion frequencies and radio flux over Cycle 23, but the pattern 
changed since the minimum of Cycle 23 with larger frequency 
shifts for a given level of radio flux during the rising phase of 
Cycle 24. The middle and right panels of Fig. [3] show respec¬ 
tively the low- and the high-frequency shifts as a function of 
the corresponding radio flux. The different behavior between the 
lower and the higher frequency modes is indicating differences 
in the structural and magnetic changes as a function of depth in 
the sub-surface layers of the Sun between Cycle 23 and Cycle 
24 (see also the discussion in Section 4 and Appendix B for 
additional figures of hysteresis at each 1). The frequency shifts 
are actually measured to be ahead of the radio flux by about 
60 days in the velocity GOLF data, while a longer delay of about 
90 days is observed with the data from the three Sun photome¬ 
ters composing the V ariability of Solar Irr adiance and Gravity 
Oscilla tions (VIRGO; 


SoHO (ISalabert et al.L 


Frohlich et al.L Il995h instrument onboard 


MU. 


Table |2] gives the differences in the frequency shifts av¬ 
eraged over I = 0, 1, and 2 modes between the maximum 
and minimum of each solar cycle for the different frequency 


ranges as illustrated in Fig. |2] The periods of minimum and 
maximum of activity were defined by selecting the data points 
included within 15% from the minimum (or the maximum) 
of the radio flux. The mean values of the frequency shifts 
over these selected periods were obtained before calculating 
the differences. The last column in Table |2] shows the ratio 
of amplitudes between Cycle 24 and Cycle 23, in the sense 
(MAX 24 -M 1 N 23 ) - (MAX 23 -M 1 N 22 ). This indicates how the 
frequency shifts respond to the weaker Cycle 24 compared to 
Cycle 23. The low-frequency modes have comparable shifts 
between Cycle 23 and Cycle 24, while the mid- and the high- 
frequency modes are about 30% smaller during Cycle 24, which 
is comparable to the decrease observed in radio flux emission 
between Cycles 23 and 24. The differences between different 
frequency ranges thus suggest that the changes associated to the 
weaker Cycle 24 occur in shallower layers of the Sun, while 
the structure and the magnetic field deeper than 1400 km must 
have remained unchanged between Cycle 23 and Cycle 24, as 
inferred from the low-frequency modes. 


3.3. At individual angular degrees 

Oscillation modes with d ifferent angular degre es I respond dif¬ 
ferently to the solar cycle (ISalabert et alilMCM in r elation to the 
spatial distribution of the s urface magnetic field (iHowe et al.L 
I 2 OO 2 I: iJimenez-Reves etML 120041) . Studying them individually 
can provide insights about the differences observed between ac¬ 
tivity cycles (see also Appendix B). 

The twelve panels of Fig. |4] show the temporal variations of 
the frequency shifts extracted from the GOLF data at each in¬ 
dividual I and calculated over the four frequency ranges. The 
individual I = 0, 1, and 2 modes are represented from top to 
bottom, while the different frequency ranges are shown from left 
to right as illustrated in Fig. [T] The 10.7-cm radio flux is also 
represented, scaled to match the rising phase and the maximum 
of Cycle 23 as previously. The corresponding gradients and the 
correlation coefficients between frequency shifts at each / and ra¬ 
dio flux are also given in Table [1] Each individual I modes show 
different responses depending on the frequency range, indicat¬ 
ing different temporal and latitudinal sensitivities to the struc¬ 
tural and magnetic sub-surface changes with solar activity. The 
gradients between each I are of the same order of magnitude 
over a given frequency range. Nevertheless, the detailed tempo¬ 
ral evolution of the individual frequency shifts is different from 
one degree to the other, and from one frequency range to the 
other, providing precious insights into the structural and mag¬ 
netic changes as a function of latitude and depth that could be 
related to the dipole and the quadrupole components of the mag¬ 
netic field, as discussed in the next section. 

4. Discussion 

The temporal variability of the solar oscillation frequencies with 
solar activity results from the combination of structural changes 
in the sub-surface layers of the Sun, associated to other changes 
resulting from the topology of the magnetic field. In fact, the in¬ 
ternal magnetic field does not act directly on the structure as the 
magnetic pressure is very low compared to the gas pressure in 
the solar interior, so it does not perturb it directly. This is the 
reason why the ID modeling of the Sun, through the Standard 
Solar Model (SSM) or the Seismic Solar Model (SeSM), do 
not contain any equations including the magnet ic field (e.g., 
iTurck-Chieze & Couvidalll201 lUBasu et al.L|2014l) . 
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Fig. 3. (Left panel) Frequency shifts, (Av„_/=o,i, 2 ), in /rHz averaged over 15 radial orders over 1800 /rHz < v < 3790 yuHz as 
a function of the corresponding 10.7-cm radio flux, F\qj, once the QBO’s signature was removed. The rising (black dots) and 
declining (blue dots) phases of Cycle 23, and the rising phase (red dots) of Cycle 24 are indicated by different colors. (Middle 
panel) Same as left panel, but for the low-frequency modes averaged between 1800 //Hz < v < 2450 //Hz. (Right panel) Same as 
the left and middle panels, but for the high-frequency modes between 3110 //Hz < v < 3790 //Hz. 

Table 2. Peak-to-peak variations of the frequency shifts, (Av„j=o,i, 2 ). averaged over the modes / = 0, 1, and 2 (in nHz) calculated 
between different activity phases of Cycles 23 and 24 over four different frequency ranges. 


Frequency range 

(MAX 23 -M 1 N 22 ) 

nHz 

(M 1 N 23 -MAX 23 ) 

nHz 

(MAX 24 -M 1 N 23 ) 

nHz 

(MaX24 -M 1 N 23 )/(MaX23 -M 1 N 22 ) 

1800 //Hz < y < 3790 //Hz 

142.9 ± 5.3 

-168.1 ±3.9 

114.1 ±4.3 

0.80 ± 0.04 

1800 //Hz < y < 2450 //Hz 

36.4 ± 7.4 

-54.3 ± 5.7 

39.9 ± 6.5 

1.09 ±0.29 

2450 //Hz<y<3110//Hz 

198.8 ±9.0 

-220.6 ± 6.4 

131.7 ±6.8 

0.66 ± 0.05 

3110//Hz<y< 3790//Hz 

455.6 ± 14.5 

-500.1 ± 10.5 

281.5 ± 10.9 

0.62 ± 0.03 


It is now well established that the variability of the acous¬ 
tic modes reveals a very small but important zone of the Sun 
located just belo w the surface called th e superadiabatic region 
- see Figure 1 of iLefebvre et al.l (l2009h illustrating the stratifi¬ 
cation of the upper laye rs of the Sun. The 3D sim ulations per¬ 
formed with STAGGER (iNordlund & Steinill998h represent one 
of the best ways to study the action of the magnetic field in this 
region as these simulations correspond to layers located between 
-2000 km and -b900 km around the surface. STAGGER properly 
describes the turbulent pressure and the motions that act below 
the solar surface, and it reproduces remarkably well the convec¬ 
tive motions at the surface. The addition of a magnetic field has 
two effects in that region of the Sun; it partly inhibits the turbu¬ 
lent pressure which cannot be properly described in a ID solar 
model, and it directly acts on the gas pressure when the magnetic 
pressure is relatively not negligible around the solar surface. The 
advantage to combine a re al 3D turbulent s tructure with a ID 
model was demonstrated bv IPiau et al.l(l2014t) . as it improves the 
prediction of the absolute values of the acoustic frequencies and 
help to interpret the observations of the sub-surface regions. The 
introduction of an accurate turbulent pressure modifies the ab¬ 
solute value of the theoretical acoustic mode frequencies. While 
the SSM or SeSM lead to theoretical frequencies greater than 
the observed ones by up to 10 //Hz, an hydrodynamical model 
reduces the differences with the seismic values within 3 //Hz be¬ 
cause these layers are better described. Eigure|5]shows this result 
compared to different ways to treat the convection in ID. 

The introduction of a magnetic field strength of 1200 G at 
the bottom of the simulation compatible with the strength of 
the toroidal field at 0.996 Rq extracted bv iBaldner et al.l ( 20091) 
from the seismic observations collected by the Michelson 


Doppler Imager (MDI; IScherrer et al.l Il995h onboard SoHO, 
reduces slightly the turbulence and leads to a small s hrink in the 
solar radius of dozens of kilometers at the surface (iPiau et al.L 
1201 ih . The I — 0 radial differences at the high-frequency 
range extracted from the GOLE observations (Eig. HJi of about 
1.3 X 10^“* between minimum and maximum of activity could be 
attributed to pure structural changes. If we follow the relation 
X “ T would correspond to an upper limit for a change 

in solar radius of about 50 km at maximum of Cyc le 23 and 
of about 30 km during Cycle 24. iBush et ^ (1201 5h extracted 
a radius change of 4 x 10“^ during the last 4 years using 
the obs ervations from the H elioseismic and Magnetic Imager 
(HMI; IScherrer et al.L l20I2h instrument onboard the Solar 
Dynamics Observatory spacecraft, which is compatible with the 
estimate extracted from the radial mode observed by the GOLE 
instrument. The / = 1 dipole and I - 2 quadrupole temporal 
variations shown in Eigure 4 at the high-frequency range, 
show reasonably well the same order of changes combined 
with a clear increase of the quasi-biennal component which 
probably translates the beating of the dipolar and qua drupolar 
components of the magnetic f i eld just below the surface (iTobiasl 

ll996tlsimoniello et al.l 1201 jtiTobias & Cattaneo[l2013h. 

The temporal variations of the medium-frequency modes 
are clearly less sensitive to rad ius changes (0. 1 5 //H z instead 
0.4 //Hz for the / = 0 modes). ILefebvre et al.l (l2009l) demon¬ 
strated that a change in radius of 2 x 10 perturbs the inter¬ 
nal structure down to 0.99 R© through opacity and equation 
of state variations that could induce pulsations of these layers 
by kappa mechanisms. The differences observed here cannot be 
then interpreted as a pure structural change when the magnetic 
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Fig. 4. Same as Fig. [1] but, from top to bottom, for the temporal variations in pHz of the frequency shifts, (Av„_;=o) , {Av„,;=i) , 
(Av„j= 2 ), at each individual angular degree / = 0, 1, and 2 (black dots). The associated 10.7-cm radio flux, F 10 . 7 , scaled to match the 
rising phase and the maximum of Cycle 23 is shown as a proxy of the solar surface activity (solid line). 


held in the sub-surface layers slightly changes. In order to bet¬ 
ter understand the situation, we also need to consider the ob¬ 
servations of the emerging solar magnetic held. Three distinct 
results are useful to that respect: 1 ) the observation of the so¬ 
lar wind collected by the Ulysses mission has shown very dif¬ 
ferent minima for Cycles 22 and 23, the minimum of Cycle 23 
having a rather opened held near the equator w hich shows the in¬ 
creasi ng importance of the toroidal component (iMcComas et al.L 
1200 ^ : 2 ) the dipolar and the radial ma gnetic fields obser ved at 
the Wilcox Solar Observatory (WSO) ilHoeksemn l2009h since 
1978; and 3) the measurements of the sunspot magnetic held 
strength which is decreasing since Cycle 22 and the measure¬ 
ment of the 1500 G threshold n eeded in order to dark sunspots 
to form (iLivingston et al.L[2012l) . The variability of the I - 1 and 
I - 2 modes is practically twice larger than the I - 0 mode vari¬ 
ability, with a more perturbed temporal evolution than the one 
observed in the 10.7-cm radio flux and the I - 0 modes. This is 
particularly true for the I = 2 modes. This is actually in agree¬ 
ment with the observations of the magnetic held performed at 
WSO which indicate that even if both dipolar and quadrupolar 
components of the external fields are decreasing since at least 
1985, the ratio between quadrupolar and dipolar components o f 
the radial held shows a constant increase (linceoglu et aLLl2015l) . 


The low-frequency range is the most interesting region to 
consider. Even if the effects are very small and that the maxi¬ 
mum of Cycle 23 is visible, we observe very perturbed signals 
for the low-degree modes mainly for the 1 = 2 modes. Table [1] 
shows low levels of correlation with the 10.7-cm radio flux since 
the maximum of Cycle 23, and Table |2] indicates that the ratio 
of the peak-to-peak variations of the frequency shifts between 
Cycle 24 and Cycle 23, (MAX 24 -M 1 N 23 ) - (MAX 23 -M 1 N 22 ), is 
larger for the low-frequency modes than for the medium- and 
high-frequency modes. Moreover, the toroidal held, responsible 
for the local partial inhibition of the turbulence and associated 
consequences, was observed to increase by about a factor 3 dur¬ 
ing the maximum of Cy cle 23 toward the sol ar interior between 
0.999 Rq and 0.996 Rq (iBaldner et al.Ll200^ . What is unclear is 
to know if the observed variability of the acoustic low-frequency 
modes comes from oscillations due to a kappa mechanism or 
from a beating between the poloidal and toroidal components of 
the magnetic held. The first assumption implies that the QBO’s 
signature should decrease as a function of time, which is not 
the case. We can then speculate about an amplification of the 
variations of the low-frequency shifts to an increase of the ra¬ 
tio between components of the magnetic held from Cycle 22 
to Cycle 24. Nevertheless, the sudden increase of the I = 1 
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shifts at low frequency is unexpected in this context and justifies 
then more investigations (see Appendix B). Also, as observed in 
sunspot number, the last solar cycles present clear hemispheric 
asymmetric^ with successive dominant north and south hemi¬ 
spheres during the maximum of activity. 

The study presented here indicates that the deeper toroidal 
magnetic field seems rather stable along Cycle 23 and 24, 
while the magnetic field at the surface is overall decreasing. 
As a consequence, a relative disconnection between acoustic 
frequencies and surface proxy, such as the 10.7-cm radio flux, is 
observed to be more and more visible in the deeper sub-surface 
layers of the Sun below 1400 km. Its physical reason is not yet 
understood. One possibility is the existence of two dynamos, a 
deep one and a shallow one. Another possibility, among others, 
is that some variability, larger than the 11-year timescale, exists 
in the tachocline and produces the emergence at the surface 
of varying flux tubes in addition to the small-scale convection 
present in the sub-surface layers. Moreover, the chao tically 
modulated solar dynamo, proposed bv lTobias et al.l (Il995h . does 
not appear to be in contradiction with the present observations. 
Any progress supposes a good knowledge of the deep magnetic 
field in the quiet and active Sun. More realistic 3D spherical 
simulations of the global dynamo coupled to superficial carte¬ 
sian 3D simulations with different topologies of the magnetic 
field could provide invaluable understanding of the internal 
dynamical Sun. 


5. Conclusions 

We analyzed 18 years of the Sun-as-a-star, space-based radial 
velocity helioseismic GOLF observations spanning the entire 
Cycle 23 and the rising phase of the weak Cycle 24. We stud¬ 
ied the temporal variations of the low-degree oscillation fre¬ 
quencies along this 18-year long dataset in order to propose 
inferences on the structural and magnetic changes in the sub¬ 
surface layers of the Sun which could explain the progression 
of Cycle 24. Indeed, the averaged kernels of the acoustic modes 
at different frequency ranges have their largest sensitivities at 
different depths. While the sensitivity of the low-degree, high- 
frequency modes (between 3110 yuHz < v < 3790 yuHz) peaks at 
0.9989 Rq (i e., 760 km) from the Sun’s surface, the low-degree, 
low-frequency modes (1800 yuHz < v < 2450 //Hz) are mostly 
sensitive to deeper-lying layers at 0.9965 Rq (i.e., 2400 km). 

We showed that while the frequency shifts associated to lay¬ 
ers above 1400 km have been following the decrease of surface 
activity during Cycle 24 compared to Cycle 23, the frequency 
shifts associated to layers below 1400 km remained similar be¬ 
tween Cycle 23 and Cycle 24 for the extrema of the Schwabe 
cycle. The study of the temporal variations of the individual 
mode oscillation frequencies observed by GOLF thus suggests 
that while Cycle 24 is clearly magnetically weaker in the shal¬ 
lower parts of the Sun, the solar Cycle 23 and 24 are magneti¬ 
cally and structurally similar below 1400 km, except the quasi¬ 
biennial oscillation for the radial and quadrupole modes which 
is observed to be slightly amplified and to possibly have a longer 
temporal period in these deeper-lying parts during the rising 
phase of Cycle 24. The dipole modes show however a specific 
and not explained behavior yet. Inversions combining low- and 
high-degree modes would provide a more precise localization of 


^ Hemispheric sunspot numbers are provided by SILSO, Royal 
Observatory of Belgium, Brussels at http://sidc.oma.be/silso/datafiles 


the evolution of the magnetic field from the surface to deeper 
parts of the Sun. 
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Fig. 5. Differences in yuHz between observed GOLF frequencies used for sound speed inversion (iTurck-Chieze & Lonesl 1201 2h 
and theoretical frequencies of t he SSM scaled by the inert ia ratio Q„ i using MLT convection (diamonds); the Canuto prescrip- 
tion of convection (trian gles) (ICanuto & Mazzitd IiL Il99 lh: t he mean thermal description (crosses); and 3D simulations from 
(iNordlund & SteinLll998l) (sauaresl. See also lPiau et alT i 20141) for the results obtained with the MDI data. 
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Appendix A: GOLF/SoHO low-degree acoustic 
oscillation frequencies 

Table A.l illustrates the content of the frequency tables made 
electronically available through this work. They contain the cen¬ 
tral frequencies of the low-degree acoustic modes of oscillations 
of the Sun extracted from 365-day sub series of 18 years of 
GOLF/SoHO observations between April 11,1996 and March 5, 
2014. As a four-time overlap of 91.25 days was used, a total of 
70 frequency tables are provided. Note then that one every four 
frequency tables contain frequencies extracted from independent 
sub series. See Section 2 for a detailed description of the method¬ 
ology applied to extract the mode frequencies. Quality criteria 
were dehned based on the htted mode parameters and their as¬ 
sociated uncertainties in order to remove outliers, such as: 1) the 
error of the mode frequency must be less than its mode width; 2) 
the signal-to-noise ratio must be larger than 1; and 3) the mode 
width must be larger than the frequency resolution. The entire 
set of tables containing the GOLF / = 0, 1, 2, and 3 p-mode fre¬ 
quencies between 1500 /iHz and 4000 /rHz, and their associated 
formal 1 cr uncertainties are available in a machine readable for¬ 
mat at CDS. The GOLF velocity time series used in this work, 
as well as the extracted frequency tables are directly available on 
the GOLF website at CEAQ and through the SoHO data archive 
at NASA and ESA. 


The GOLF website at CEA is also mirrored through the SPACEINN 
website http://www.spaceinn.eu 


Appendix B: Frequency shifts at individual angular 
degree 

Eigure lB/Tl shows the temporal variations of the frequency shifts, 
(Av„_;=o) , (Av„,/=i) , (Av„,/= 2 ), at each individual angular de¬ 
gree 1. The shifts were calculated over four frequency ranges: 
1800 //Hz < V < 3790 //Hz; 1800 //Hz < v < 2450 //Hz 

(the low-frequency range); 2450 //Hz < v < 3110 //Hz (the 
mid-frequency range); and 3110 //Hz < v < 3790 //Hz 

(the high-frequency range). The QBO’s signature was removed 
by applying a prope r smoothing of 2.5 years as prescribed in 
lEletcher et al.l (1201 Ol) for the four frequency ranges and for the 
10.7-cm radio flux. Note that in order to avoid misinterpreta¬ 
tion due to possible border effects introduced by the smooth¬ 
ing, we do not consider the extreme points corresponding to 
the hrst three and the last three points of the smoothed data. 
Each frequency range is sensitive to differen t depths in t he sub - 
surface layers of the Sun, as calculated by iBasu et aTI (l2012l) . 
The averaged kernels of the mid- and the high-frequency ranges 
have their largest sensitivities at 0.9981 Rq (i.e., 1300 km) and 
0.9989 Rq (i.e., 760 km) respectively, while the low-frequency 
modes peak deeper in the solar interior at 0.9965 Rq (i.e., 
2400 km). See Section 3 for results about the frequency shifts 
averaged over the modes / = 0, 1, and 2, (Av„j=o,i, 2 )- 

While the temporal variations of the frequency shifts at / = 0, 
1, and 2 show a good agreement with the 10.7-cm radio flux for 
the mid- and high-frequency modes, clear differences between 
the two quantities can be seen for the low-frequency modes. A 
residual after removing the QBO is still present in the I - 0 and 
/ = 2 modes. The longer period observed in the low-frequency 
modes starting around 2006 and mentioned in Section 3 appears 
in the I = 0 and I = 2 modes only. In the other hand, the low- 
frequency I = 1 has a very different behavior. After a long period 
between 2004 and 2010 of practically no changes one the fre- 
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Table A.l. Acoustic oscillation frequencies of the / = 0, 1, 2, and 3 modes extracted from the 365-day GOLF/SoHO spectrum 
calculated between 11-04-1996 and 10-04-1997. 


n 

/ = 0 

/ = 1 

1 = 2 

1 = 3 

9 

+ 

± 

1535.872 ± 0.034 

± 

10 

1548.414 ± 0.041 

1612.734 ± 0.028 

1674.568 ± 0.057 

± 

11 

1686.538 ± 0.043 

1749.297 ± 0.040 

1810.304 ± 0.052 

+ 

12 

1822.197 ±0.043 

1885.106 ±0.037 

1945.823 ± 0.045 

+ 

13 

1957.469 ± 0.051 

2020.840 ± 0.033 

2082.058 ± 0.049 

2137.784 ± 0.064 

14 

2093.510 ±0.038 

2156.811 ±0.053 

2217.819 ± 0.083 

2273.370 ± 0.075 

15 

2228.819 ± 0.049 

2291.996 ±0.064 

2352.285 ± 0.074 

2407.785 ±0.122 

16 

2362.799 ± 0.065 

2425.473 ± 0.060 

2485.921 ± 0.083 

2541.352 ±0.126 

17 

2496.222 ± 0.060 

2559.231 ±0.059 

2619.644 ± 0.070 

2676.191 ±0.098 

18 

2629.906 ± 0.055 

2693.303 ± 0.052 

2754.474 ± 0.061 

2811.351 ±0.087 

19 

2764.146 ±0.052 

2828.103 ±0.056 

2889.557 ± 0.053 

2946.983 ± 0.058 

20 

2898.947 ± 0.055 

2963.273 ± 0.054 

3024.664 ± 0.064 

3082.219 ± 0.093 

21 

3033.772 ± 0.054 

3098.152 ±0.063 

3159.856 ±0.066 

3217.738 ±0.088 

22 

3168.619 ±0.055 

3233.168 ±0.065 

3295.107 ±0.090 

3353.672 ±0.164 

23 

3303.376 ±0.071 

3368.612 ±0.081 

3430.871 ±0.109 

3489.737 ± 0.361 

24 

3438.884 ± 0.086 

3503.887 ±0.105 

3566.714 ± 0.200 

3626.576 ± 0.466 

25 

3574.584 ±0.170 

3640.344 ±0.135 

3703.153 ±0.348 

± 

26 

3710.794 ± 0.235 

3776.536 ± 0.208 

3839.292 ± 0.638 

± 

27 

3847.099 ± 0.334 

3913.674 ± 0.267 

+ 

± 


quency shifts, a sharp increase is measured starting from 2011, 
which is actually not present in the other angular degrees. 

Figure ED shows the same frequency shifts, {Av„j=o) , 
(Av„,/=i) , (Av„,/= 2 ), at each individual angular degree I, but as 
a function of the corresponding 10.7-cm radio flux, T^in ?. Bot h 
quantities were properly smoothed as in iFletcher et al.l (l2010h . 
The shifts obtained over three frequency ranges are represented; 
1800 pHz < V < 3790 jMz; 1800 /Mz < v < 2450 fiUz (the 
low-frequency range); and c) 3110 yuHz < v < 3790 pHz (the 
high-frequency range). Different colors are used for three dif¬ 
ferent phases of solar activity: a) the rising phase of Cycle 23 
going from April 1996 to October 200; b) the declining phase of 
Cycle 23 going from October 2001 to January 2009; and c) the 
rising phase of Cycle 24 going from January 2009 to April 2014 
(see Section 3 for more details). Again, the low-frequency modes 
show a very different behavior compared to higher-frequency 
modes. Moreover, a longer period than the quasi-biennial oscil¬ 
lation is observed for the radial and quadrupole modes during 
the rise of Cycle 24, while the dipole modes show a very distinct 
behavior. 
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Fig. B.l. Temporal variations in fiHz of the frequency shifts, (Av„_/=o) , {Av„_;=i) , (Av„_;= 2 ), at each individual angular degree I 
(from top to bottom), and calculated for four different frequency ranges, the QBO’s signature being removed (black dots). From 
left to right, the frequency ranges are the following: a) 1800 //Hz < v < 3790 //Hz; b) 1800 //Hz < v < 2450 //Hz; c) 2450 //Hz 
< V < 3110 //Hz; and d) 3110 //Hz < v < 3790 //Hz. The 10.7-cm radio flux, F 10 . 7 , smoothed in the same way and scaled to match 
the rising phase and maximum of Cycle 23, is shown as a proxy of the solar surface activity (solid line). 
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Fig. B.2. Frequency shifts, (Av„,/=o) . (Av„_;=i) , {Av„j= 2 }, at each individual angular degree / = 0, 1, and 2 (from left to right), 
as a function of the corresponding 10.7-cm radio flux, Fio.?, and calculated for three different frequency ranges once the QBO’s 
signature was removed. From top to bottom, the frequency ranges are the following: a) 1800 yuHz < v < 3790 //Hz; b) 1800 //Hz 
< V < 2450 //Hz (the low-frequency range); and c) 3110 //Hz < v < 3790 //Hz (the high-frequency range). The rising (black dots) 
and declining (blue dots) phases of Cycle 23, and the rising phase (red dots) of Cycle 24 are indicated by different colors. 
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